Nature Methods

Subcellular mass spectrometry imaging (MSI) has required specialized instrumentation and approaches. Additionally, data analysis has required expensive proprietary software, or code-based open-source tools. These barriers render spatial metabolomics less accessible than peer spatial omics techniques (transcriptomics, proteomics). Here we demonstrate subcellular MSI on standard commercial hardware, analyzed using code-free, open-source software. An inexpensive, fully reversible modification of a stock DESI source takes the platform past its typical 5-10 m pixels to routine 2x2 m pixels, and 1x1 m in proof-of-principle experiments. At a tissue-glass boundary imaged at 2x2 m, tissue-specific ions showed no measurable spillover onto the adjacent glass within the {+/-}4.4 m uncertainty of the edge position. To make these data interpretable, we extend MSI.EAGLE, an open-source, vendor-agnostic application with two new methods: i) Spatial UMAP preserves tissue architecture alongside subcellular detail, where conventional clustering breaks apart at high resolution; ii) Histology Fit performs landmark-free co-registration, and recovers synthetic misalignments of up to 8 m to within one pixel of a reference registration. Together, these enable previously high-burden spatial metabolomics experiments. In mouse brain tissue, we demonstrate genuine subcellular signal by resolving mitochondrial cardiolipin into perinuclear puncta, quantifying cytoplasmic enrichment across 6,073 cells. Cell identities are assigned to MSI pixels using histology-predicted transcriptomics, linking cell provenance to metabolism in 2D space directly. Finally, we recover cell-type-specific metabolomes consistent with the transcriptomic identities. This workflow allows cellular and subcellular spatial metabolomics linked to cell provenance, with only stock MSI hardware, open-source GUI-based software, and common histology techniques, putting it within reach of many researchers and laboratories.

Motion artifacts are almost unavoidable in functional time-lapse and structural volumetric multiphoton microscopy. They arise from respiration, heartbeat, locomotion, awake behavior, instrument heating, mechanical vibration, and slow drift, while the recorded signal is often photon-limited, blurred by scattering, and biologically time varying. Consequently, motion correction is frequently an essential prerequisite for quantitative bioimage analysis rather than a merely cosmetic preprocessing operation. Edge- and landmark-centric registration strategies are often poorly matched to these data because useful structures may be sparse, diffuse, out-of-focus, or changing in fluorescence intensity. We present ZenReg, an open-source Python platform that formulates common 2D+t, 3D, and 3D+t microscopy registration tasks as a modular family of geometry-preserving alignment problems. ZenReg combines Fourier phase correlation, intensity-based StackReg-style alignment, NoRMCorre-style piecewise translation fields, projection-based rotation estimates, dense SimpleITK-based six-degree-of-freedom volume registration, and sparse point-based rigid-volume registration within one canonical microscopy stack model. The platform uses OMIO to normalize heterogeneous microscope files and to preserve metadata, while optional disk-backed Zarr arrays support chunked, memory-efficient processing of image stacks that exceed available memory or reside on remote storage. ZenReg writes registered images together with shift tables, correlation metrics, summary plots, and machine-readable settings. In synthetic benchmarks with known ground truth, ZenReg recovered global 2D and 3D translations with subpixel accuracy across moderate noise and drift regimes, while high-noise and large-drift tests separated the backend behavior: FFT-based methods failed abruptly once image information or shared support became insufficient, StackReg degraded more gradually under severe noise, and piecewise NoRMCorre improved spatially varying local-motion correction where a single global transform was inadequate. Parallel execution reduced runtime for large time series, and ZenReg provided practical full-volume rigid correction for dense and puncta-rich 3D+t stacks. By coupling a modular, extensible backend architecture to transparent sidecar outputs, ZenReg makes motion correction easier to extend, inspect, share, reproduce, and reuse as part of scientific image analysis.

Light microscopy is routinely used to explore the cellular and molecular structure of tissues, but the scale and complexity of data remains a bottleneck for discovery. We introduce 3D-MAESTRO (3D-Microscopy Automation and Execution with Scalable Tools, Rendering, and Orchestration): an automated image processing workflow for large-scale microscopy data, built for scalable execution across cloud and local computing environments. 3D-MAESTRO orchestrates denoising, stitching of image tiles, atlas registration, and segmentation. Its modular architecture permits the integration and benchmarking of new packages, ensuring that performance evolves as more efficient or accurate algorithms emerge. We introduce a new 3D image template for automated registration of mouse brains cleared with aqueous reagents and an efficient method for detection of fluorescent cells. We apply 3D-MAESTRO to lightsheet images of whole mouse brains in the context of diverse anatomical and functional experiments, illustrating high-throughput and reproducible mapping of microscopic structures across the brain.

We developed EASI-PASS, a reliable, high-throughput method for estimating the molecular identity of functionally characterized cells by merging live imaging with subsequent fixed-tissue imaging using conventional microscopes. Our method matches the shapes and locations of thousands of densely imaged cells between large (>1 mm2) functional images and a thick, expanded, and cleared EASI-FISH tissue volume to assess gene expression. This approach is more efficient than alignment to thin sections and recovers the molecular identity of ~78% of cells. In acute brain slice imaging from the mouse parabrachial nucleus during optogenetic stimulation of long-range spinal inputs, we observed fine-scale specificity in the molecular identity of spinorecipient neurons. In the awake mouse visual cortex, we observed distinct arousal modulation and spatial falloff in correlations within and across interneuron classes. Thus, EASI-PASS provides reliable and efficient alignment of cellular activity with molecular identity.

Understanding distributed biological systems, particularly neural circuits, requires simultaneous cellular-resolution imaging across millimeter-scale fields of view (FOV). Existing miniature microscopes remain fundamentally constrained by trade-offs among FOV, spatial resolution, and optical complexity, limiting their ability to bridge cellular microscopy with cortex-scale imaging. Here we introduce distributed computational optics, a framework that distributes image formation across coordinated optical modules and computationally integrates their measurements into a unified image. We realize this framework in Bio-CM2, a computational miniature mesoscope that partitions the imaging field across four optical modules while converging their measurements onto a common image sensor. This architecture overcomes the aberration-scaling limitations of conventional miniature optics while avoiding the hardware complexity of multi-camera systems and the contrast degradation associated with optical multiplexing. Bio-CM2 achieves a 7.5 x 10 mm2 FOV while enabling cellular-resolution in vivo imaging at video rates. We demonstrate its utility through two complementary imaging modalities in head-fixed mice: cortex-wide functional vascular imaging, enabling simultaneous quantification of pial arteriole vasomotion and mesoscale hemodynamic functional connectivity, and cellular-resolution calcium imaging, resolving the activity of over 3,000 neurons together with mesoscale neuronal functional connectivity. We further demonstrate the versatility of the platform through cellular-resolution imaging of entire coronal mouse brain sections, population-scale imaging of freely behaving Caenorhabditis elegans, and odor-evoked calcium imaging of the main olfactory bulb in head-fixed mice, highlighting its broad applicability across diverse biological systems and imaging modalities. By overcoming the conventional trade-off between FOV and spatial resolution in a compact miniature platform, Bio-CM2 establishes distributed computational optics as a scalable framework for multiscale biological imaging.

Biological image analysis is full of discrete decisions: whether an object is present, which of several overlapping detections is real, whether two detections match across time, and whether a cell divides. Standard pipelines resolve them locally with non-max suppression, thresholding, or greedy linking, committing before all image and temporal evidence is in. We recast such problems as discrete inverse rendering: candidate renderings are generated then jointly selected to reconstruct the movie subject to temporal and biological constraints, solved to certified optimality with a modern integer-programming solver. The same formulation covers suppression of overlapping detections, selection of a structure as a path, and event-structured tracking with birth, death, and division. Applied to C. elegans splines, sperm flagella, and dividing cells, the method matches specialised state-of-the-art pipelines across three imaging modalities on a single objective, with the largest gains where per-frame segmentation is unreliable (on a low-signal fluorescence movie of Huh7 hepatoma cells, detection F1 doubles from 0.31 to 0.58). HighlightsO_LIInteger-programming framework for suppression, path, and lineage selection C_LIO_LIOne objective, one solver: worm splines, sperm flagella, and dividing cells C_LIO_LIReconstruction-based scoring matches specialised pipelines across three modalities C_LIO_LICertified-optimal solutions in seconds to minutes on standard benchmark movies C_LI In BriefZdyb and Kirkegaard recast several biological image-analysis problems--non-max suppression, structural path extraction, and event-structured cell tracking--as one discrete inverse-rendering problem solved by an integer-programming solver, matching specialised trained pipelines across three distinct imaging modalities within a single objective.

Calcium transients encode cellular and neuronal activity across timescales ranging from milliseconds to hours, yet linking these transient signals to downstream molecular states remains a major challenge. We recently introduced Caprola, a calcium-dependent protein labeling tool that converts calcium transients into permanent fluorescent marks for later analysis. In this way, Caprola enables tracking of neuronal activities in animal models as well as retrospective identification of labeled cells for isolation and transcriptomic analysis. However, the relatively slow labeling kinetics of Caprola required high concentrations of fluorophore probe and relatively long labeling times, which limits its sensitivity and applicability, in particular in vivo. To address this limitation, we generated Caprola variants with up to 29-fold faster labeling rates than their predecessor. We demonstrate that our new Caprola variants record calcium transients in cells and in zebrafish larval brains under conditions where previous Caprola variants did not show labeling. We further expand the applicability of Caprola to activity-dependent marking of postsynaptic compartments, opening new avenues for coupling functional activity histories with downstream molecular and transcriptomic analyses.

Structural characterization of macromolecular assemblies within intact cells remains one of the central challenges in cell biology. While cryo-electron tomography provides unparalleled structural information, its applicability is limited by sample thickness, imaging throughput, and accessibility. Fluorescence microscopy offers molecular specificity and compatibility with intact biological specimens but has so far lacked the spatial resolution required to visualize cellular ultrastructure. Here we introduce Mega-expansion microscopy (Mega-ExM), a fluorescence imaging approach that enables structural visualization of whole cells using conventional confocal microscopes. Mega-ExM combines iterative hydrogel expansion with whole-proteome NHS-dye labeling and post-expansion immunostaining to achieve tunable expansion factors of up to [~]1,500-fold while preserving ultrastructure. At expansion factors of 40-260x, Mega-ExM resolves centrioles, mitochondrial cristae, protein-dense domains within mitochondrial cristae consistent with respiratory-chain supercomplexes, the synaptonemal complex, and nuclear pore complexes (NPCs) with high fidelity. Particle averaging of [~]200x expanded NPCs yields reconstructions with a structural resolution of [~]35 [A], approaching what cryo-electron tomography has achieved for selected protein assemblies. By combining molecular specificity, large imaging volumes, and nanoscale structural resolution on conventional fluorescence microscopes, Mega-ExM establishes a broadly accessible platform for in situ structural biology.

Single-cell spatial transcriptomics can connect molecular cell states with tissue morphology, but this promise depends on accurate registration to histopathology. In serial sections, however, tissue borders often differ because of sectioning artifacts, staining variability, and field-of-view acquisition, limiting conventional area-based registration. We developed AnchorR, an expert-guided workflow for coarse-grained alignment of hematoxylin and eosin (H&E) images with CosMx Spatial Molecular Imaging data. Bioinformaticians first define and color-code cell types in Seurat, and pathologists then identify corresponding internal landmarks using QuPath overlays. AnchorR combines these paired landmarks to estimate affine transformations, quantify residual error, and support visual quality control and anchor refinement. Using six oral pre-cancerous tissue sections, we identified 60 cross-modal landmarks. Fitting each section independently reduced mean landmark error from 121.5 {micro}m with a single whole-slide transformation to 14.6 {micro}m. Cross-validation further showed that increasing the number of anchors improved robustness, with nine-anchor fits achieving approximately 20 {micro}m error, or about one cell diameter. AnchorR is designed to complement automated computer-vision methods by providing reliable tissue-level alignment when border mismatch makes global registration difficult. By creating a shared workspace for pathologists and bioinformaticians, it operationalizes an expert-in-the-loop approach and makes feature-based multimodal registration accessible without specialized computer-vision expertise or high-performance computing.

Recent efforts have demonstrated that systematic analysis of the millions of metagenomes in the Sequence Read Archive (SRA) enables groundbreaking biological discoveries. However, SRA-scale comparison of whole metagenomes has remained elusive. Here, we introduce the most comprehensive resource to date for experiment-to-experiment analysis of all SRA metagenomes, comprising FracMinHash sketches (470GB), hypergen sketches (19GB), and an all-pairs Jaccard similarity matrix (90GB), which we release online at https://github.com/RolandFaure/SRA-metagenome-sketches-databases. These resources enable a wide range of analyses, which we demonstrate with several applications. First, we use the FracMinHash sketches to uncover temporal trends, finding that the amount of newly discovered microbial sequence continues to grow, rather than saturate or shrink. Second, we use the hypergen sketches, which enable nearest-neighbor search, to identify co-assembly candidates and thereby improve assembly quality. Third, we use the similarity matrix to analyze all SRA metagenomes, identifying duplicated submissions and over 50K candidates for metadata correction. Together, our resources unlock systematic analysis of metagenomic data at unprecedented scale.

Volumetric microscopy produces datasets too large to share as files and too rich to convey as still images, yet exploring them interactively means installing desktop software or building web infrastructure. Here we present Luxar, an open-source Python framework that compiles multidimensional data into compact archives that any modern browser renders from static file hosting. Fluorescence volumes are represented as mixtures of anisotropic Gaussians (Gaussian splats): the image becomes geometry that the GPU draws directly, with no voxel grid to rebuild. Splats share one scene with points, trajectories and surfaces. Blind-spot cross-validation sets the number of Gaussians without noise-free references. Across 17 volumes from 12 datasets and four modalities, this yields a median 99-fold (6- to 340-fold) size reduction at 26 to 67 dB peak signal-to-noise ratio (PSNR), in minutes per volume on one GPU. We demonstrate Luxar on a 500-timepoint embryogenesis recording, super-resolution localisations, single-cell atlases and lineage tracking, each explorable through a shareable link.

Focal adhesions (FAs) and actin stress fibers (SFs) are mechanically coupled but often analyzed separately. We present FilAdh, a unified framework for FA instance segmentation, SF instance reconstruction, morphometry, and FA-anchorage-based SF classification. Across 57 held-out HeLa, SUIT-2, and MRC-5 cell line images, FilAdh achieved the highest FA Dice, intersection over union (IoU), and recall and the highest SF clDice and skeleton precision among benchmarked methods; paired t tests with Holm correction indicated significant gains over all compared methods for FA Dice, IoU, and recall, and over classical methods for SF centerline metrics, with no significant differences detected relative to the deep-learning comparators. Conditional on expert-curated FA and SF masks, SFs were assigned to dorsal, ventral, and peripheral classes with 81.5%-87.3% accuracy. Controlled Poisson-Gaussian noise, Gaussian defocus, and point-spread-function blur revealed distinct, structure-specific errors in overlap and object-level measurements. By linking segmentation and reconstruction to count, area, length, apparent width, and structural class, FilAdh offers an interpretable view of the FA-actin structural system and makes degradation-dependent measurement biases explicit. BIGGER PICTURECells sense and respond to their physical environment through focal adhesions and actin stress fibers. The two form a mechanically coupled system, yet they are usually analyzed separately, with workflows re-tuned for each new set of images. FilAdh analyzes both within a single framework, tailoring segmentation, reconstruction, and classification to each structures geometry: compact adhesion patches versus thin, extended fiber networks. Testing on expert-annotated references and under controlled degradations that mimic low signal, focus errors, and optical blur shows that these common imaging imperfections distort adhesion and fiber measurements in different ways. The same design principle--matching the analysis to each structures shape and judging methods by how errors propagate into the final biological measurements--may benefit other bioimage-analysis problems. Released code, trained weights, and fixed data splits make all comparisons reproducible and extensible. HIGHLIGHTSO_LIFilAdh jointly quantifies focal adhesions and actin stress fibers C_LIO_LIFilAdh achieves the highest benchmarked FA overlap and comparable SF topology C_LIO_LIMatched synthetic degradations reveal geometry-specific measurement errors C_LIO_LIFA-SF spatial context enables dorsal, ventral, and peripheral SF classification C_LI

Understanding how cells establish spatial organization within tissues is a fundamental question in life sciences. While modern three-dimensional fluorescence microscopy captures large-volume tissue architecture, extracting quantitative cellular insights from complex volumetric datasets remains a major barrier. Here, we introduce FOCUS-3D, a robust, broadly generalizable volumetric cell segmentation framework built on a large, diverse manually annotated cell resource and advanced AI designs. Integrating volumetric representation learning, multi-scale feature extraction, and query-based mask prediction, FOCUS-3D achieves state-of-the-art performance across diverse species, tissues, fluorescent reporters and imaging modalities. During zebrafish (Danio rerio) development, FOCUS-3D uncovers three successive phases of notochord morphogenesis. We disentangle early motility-driven rearrangements from later cell shape remodeling and tissue repacking, and further link these morphological states to spatial and developmental transcriptional programs across independent datasets.

Super-resolution methods, such as single molecule localisation microscopy (SMLM), can determine the distribution of specific proteins at the nanometre scale. However, three-dimensional (3D) SMLM is typically limited to thin sections and requires prohibitively long acquisition times to image complete cells. Here we introduce a 3D SMLM imaging modality, Intrinsic3D, that combines the speed and optical sectioning of oblique plane microscopy (OPM) with the nanometre resolution of SMLM to determine the position of single molecule emitters throughout complete volumes more than an order-of-magnitude faster over regular DNA-PAINT. By scanning the oblique light sheet across the volume and combining localization information from two or more adjacent planes, an emitter's x, y and z position is determined, eliminating the need for post-processing steps to stitch multiple image stacks together. We demonstrate resolution and speed by imaging nuclear pore complexes in U2OS cells and Drosophila egg chambers. As a proof-of-concept, we map the 3D distribution of the cell-adhesion protein DE-cadherin in the Drosophila follicular epithelium and show that spot adherens junctions are enriched near tricellular junctions and contain more E-cadherin than the apical adherens junctions, suggesting that spot junctions play a wider role in cell-cell adhesion than previously appreciated. Finally, we examine the organization of Gliotactin, a tricellular septate junction protein relative to E-cadherin using exchange-PAINT. These investigations demonstrate that Intrinsic3D delivers rapid 3D protein mapping in whole cells. Combining Intrinsic3D with existing multiplex and resolution enhancing methods (e.g. exchange-PAINT and RESI) will enable multicolour, 3D in situ spatial proteomics in tissue samples and bridge the gap between SMLM and systems-level cell biology.

Interactive exploration and sharing of hybridization-based spatial transcriptomics datasets remain limited by large file sizes, specialized software, and computationally intensive preprocessing, restricting rapid access to publicly available resources such as the Brain Image Library. Here we present MERFISHEYES (https://merfisheyes.com), a browser-based platform for interactive visualization of single-cell and single-molecule spatial transcriptomics data across commercial and emerging high-throughput platforms. MERFISHEYES enables marker gene exploration, cross-dataset and cross-species comparison, visualization of subcellular RNA localization, and collaborative sharing through browser-accessible links generated directly from drag-and-drop datasets. We use MERFISHEYES to enable browser-based visualization of all 158 spatial transcriptomics datasets in the Brain Image Library and demonstrate its capabilities in handling high-throughput 3D data with genome-scale imputation of gene expression and DNA accessibility. MERFISHEYES provides a scalable framework for rapid exploration, sharing, and collaborative analysis of spatial transcriptomics data.

Quantitative cell colocalization in fluorescence microscopy often depends on ad hoc combinations of image loading, segmentation, region selection, manual inspection, and spreadsheet post-processing. Such workflows are difficult to transfer across projects and often obscure how intermediate results were produced. We present CellColoc, an open-source Python workflow pipeline for segmentation-based cell colocalization, single-channel segmentation, and cell feature extraction in 2D and 3D microscopy images. CellColoc provides a modular workflow layer that integrates existing segmentation backends, including Cellpose and threshold-based methods, into reusable, script-driven analyses. The package supports channel-wise backend selection, interactive or reusable regions of interest, optional third-channel occupancy and cell-positivity analysis, z-cropping and z-projection, cached post hoc refinement of Cellpose thresholds, and reanalysis after manual mask editing. Analyses are executed from concise user scripts while reusable functionality is kept in a core package. Intermediate artifacts such as ROI masks, per-channel label masks, positive-cell masks, and structured result tables are written to a standardized results directory, promoting transparent inspection, reproducibility, and FAIR-aligned reuse. Public example datasets, a synthetic benchmark, and archived software releases accompany the package. By separating reusable analysis logic from project-specific configuration, CellColoc offers an extensible foundation for community-driven microscopy workflows that need transparent per-cell overlap classification, morphology readouts, and reusable batch analysis.

Deep learning has enabled single-molecule localization microscopy (SMLM) at high emitter densities, but only for single channel systems. Here we present DECODE-Plex, a deep-learning-based framework to localize dense single molecules with overlapping point spread functions simultaneously in multiple channels. We showcase DECODE-Plex on experimental ultra-high density dual-color and 3D live-cell data. Packaged for ease of use, it will enable many groups to improve imaging speed and quality of multi-channel SMLM.

Quantitative fluorescence microscopy requires photon-efficiency, speed and accurate intensity and lifetime measurements. Time-correlated single-photon counting (TCSPC) simultaneously captures intensity and lifetime, but photon pile-up distorts both signals at high count rates, preventing fast acquisitions. Existing corrections discard photons, distort intensity, or require specialized detectors. Here we introduce PURELIGHT, an integrated hardware and software framework that simultaneously recovers undistorted intensities and lifetimes at count rates far beyond conventional pile-up limits. PURELIGHT works with hybrid photodetectors, silicon photomultipliers and photomultiplier tubes while retaining over three times more photons than alternative approaches. Using two-photon imaging, we showcase PURELIGHT's superior accuracy and spatial contrast, demonstrating video-rate subcellular lifetime imaging in awake mice, a unique lifetime-calibrated ratiometric modality and crosstalk-free temporal multiplexing. By removing the limits that have confined TCSPC to low-signal applications, PURELIGHT promotes the adoption of quantitative, photon-efficient microscopy across the life sciences.

Volumetric calcium imaging is increasingly used to capture larger neuronal populations at higher throughput, but high-speed axial sampling can compromise single-neuron identity. Here we identify cross-plane identity duplication as a structured error in volumetric imaging: anisotropic axial blurring and plane-wise functional segmentation can repeatedly detect the same neuron across adjacent planes, creating duplicate functional nodes that inflate neuronal counts and distort network phenotypes. We developed Comprehensive Label-Guided (CLG) volumetric imaging, a structural-functional calibration framework that uses nuclear labels as stable three-dimensional identity anchors for calcium signals. CLG combines nuclear labeling, deep-learning-based 3D segmentation, anatomical registration and identity-guided trace reassignment. In larval zebrafish whole-brain recordings, CLG resolved ~30,000 redundant detections and reduced estimated neuronal counts by 37-46%. In mouse visual cortex, CLG consolidated ~40% of putative duplicates and recovered over 2,000 active neurons missed by calcium-only analysis. Across baseline and perturbed conditions, calibration stabilized graph-derived measurements of hub organization, long-range correlations and network resilience. CLG therefore defines an anatomy-constrained identity-calibration layer for reliable single-neuron-resolved volumetric imaging.

Long-read RNA sequencing directly resolves the full structures of RNA transcripts. Advances in throughput now enable the generation of deeply sequenced cohorts of hundreds of samples, making joint transcript discovery across large datasets possible. However, existing transcript identification methods were designed for small datasets, which limits their applicability at this scale. Here, we present Isocall, a scalable and deterministic computational method for jointly calling transcripts from multiple PacBio long-read RNA sequencing samples. Isocall converts aligned full-length non-concatemer reads into compact per-sample transcript profiles, merges these profiles across samples, and jointly identifies known and novel transcripts supported by reads in the analyzed dataset. Filtering is tunable: presets provide coarse control and individual parameters, including relative abundance and internal priming thresholds, provide fine control. Isocall demonstrated high precision in our accuracy benchmarks, including the WTC11 SIRV spike-in controls, for which Isocall reported 0-2 false-positive transcripts per sample across the three SIRV mixes at default settings. To demonstrate scalability, we applied Isocall to 206 samples, totalling 3.5 billion raw reads, from the Human Pangenome Reference Consortium. After parallelized pbmm2 alignment and Isocall profile, the call step performed joint calling across the entire dataset in 25 minutes, using 1.3 GB of peak memory and 8 threads. Finally, in Genome in a Bottle samples with matched SNP genotypes, splice-site polymorphisms provide an additional measure of call accuracy: Isocall recovered 337 polymorphic splice sites, including a de novo donor site in BTN3A1 that corresponds to a complete isoform switch on the mutant allele.